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Abstract. Numerical solutions of 2-D steady incompressible flow inside a triangular cavity are presented. 
For the purpose of comparing our results with several different triangular cavity studies with different tri- 
angle geometries, a general triangle mapped onto a computational domain is considered. The Navier-Stokes 
equations in general curvilinear coordinates in streamfunction and vorticity formulation are numerically 
solved. Using a very fine grid mesh, the triangular cavity flow is solved for high Reynolds numbers. The 
results are compared with the numerical solutions found in the literature and also with analytical solutions 
as well. Detailed results are presented. 
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1 INTRODUCTION 

Flows inside closed geometries have always been the fo- 
cus of attention of Computational Fluid Dynamics (CFD) 
studies. As an example, the square driven cavity flow can 
be counted. Most probably the interest to the driven cav- 
ity flow is due to the simplicity of the geometry and also 
the simplicity of the application of the boundary condi- 
tions. Despite its simplicity, the driven cavity flow shows 
interesting flow features as the Reynolds number increases. 

In the literature it is possible to find many numerical 
studies on the driven cavity flow. Among numerous pa- 
pers found in the literature, Erturk et al. [6], Erturk and 
Gokcol [7], Barragy and Carey [1], Botella and Peyret [4], 
Rubin and Khosla [18], Benjamin and Denny [3], Ghia et 
al. [9] and Li et al. [12] are examples of numerical studies 
on the driven cavity flow. Burggraf [5] applied the Batch- 
elor's [2] model to driven cavity flow and obtained the 
theoretical core vorticity value at infinite Reynolds num- 
ber. Erturk et al. [6] introduced a new efficient numerical 
method and with this method, using a fine grid mesh, they 
solved the driven cavity flow at very high Reynolds num- 
bers up to i?e=21000. They showed that at high Reynolds 
number the vorticity at the center of the primary vortex 
asymptotically approaches to a value close to the theoreti- 
cal value obtained by Burggraf [5]. The difference between 
the numerical values of Erturk et al. [6] and the theoret- 
ical value can be attributed to the thin boundary layer 
developed along the solid walls and the counter rotating 
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vortices appear at the corners as the Reynolds number 
increases. 

The 2-D steady incompressible flow inside a triangular 
driven cavity is also an interesting subject like the square 
driven cavity flow. This flow was studied numerically by 
McQuain et al. [14], Ribbens et al. [17], Jyotsna and Vanka 
[11], Li and Tang [13] and Gaskell et al. [8]. After a brief 
literature survey, to the best of our knowledge, these are 
the only numerical studies found in the literature on the 
steady incompressible flow inside a triangular driven cav- 
ity. The results of these studies, however, show some dis- 
crepancies. This constitutes the main motivation of this 
study. 

Apart from numerical studies, Moffatt [15] studied the 
triangular cavity flow analytically in the Stokes regime. 
Moffatt [15] showed that the intensities of eddies and the 
distance of eddy centers from the corner, follow a geomet- 
ric sequence. Moffatt's [15] results will provide a mathe- 
matical check on our numerical results, as it was also used 
in Jyotsna and Vanka's [11] and Gaskell et al.'s [8] study. 

McQuain et al. [14] applied the Batchelor's [2] mean 
square law to triangular cavity flow and analytically ob- 
tained the inviscid core vorticity for infinite Reynolds num- 
ber. In this study we will discuss the predictions of the 
mean square law for the triangular cavity flow and its ap- 
plicability for the triangular cavity flow. 

Erturk et al. [6] stated that for square driven cavity 
flow fine grid mesh is necessary for the reliability of the 
solution when high Reynolds numbers are considered. In 
order to obtain high accurate numerical solutions, follow- 
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ing Erturk et al. [6], we will use a fine grid mesh with 
(512x512)/2 grid points. 

The aim of this study is then, to solve the 2-D steady 
incompressible triangular cavity flow using a fine grid mesh 
and present accurate solutions. We will compare our re- 
sults in detail with the numerical studies found in the 
literature. We will also compare our Stokes regime solu- 
tions with analytical results of Moffatt [15] and our high 
Reynolds number solutions with Batchelor's [2] model which 
consists of an inviscid core with uniform vorticity coupled 
to boundary velocities at the solid wall. Detailed results 
will be presented. 



1 /8 2 uj 8 2 uj\ dip 8lu dip dtv 
Re V dx 2 dy 2 J dy dx dx dy 

where ip is streamfunction and ui is vorticity, x and y are 
the Cartesian coordinates and Re is the Reynolds number. 
We note that these equations are non dimensional and a 
length scale of (h — c y )/3 and a velocity scale of U, ie. the 
velocity of the lid, are used to non-dimensionalize the pa- 
rameters and the Reynolds number is defined accordingly. 

Using the chain rule, the governing N-S equations in 
general curvilinear coordinates are the following 



2 FORMULATION AND NUMERICAL 
PROCEDURE 

Different studies found in the literature concerning flows 
in triangular cavities use different triangle geometries. For 
the purpose of being able to compare our results with 
all these studies, we have considered a general triangle 
geometry. A general triangle with corner points A, B and 
C where A and B have the same y position, as shown in 
Figure 1, can be mapped onto a computational domain, 
also as shown in Figure 1, using the following relation. 



C X H (l X Cy 
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Also the inverse transformation can be done using the 
following relation 

y = (c y - h)t] + h 



x=(b x - a x )£ - (a x - c x )rj + , 



(2) 

From these relations, we can calculate the transforma- 
tion metrics as the following 
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The 2-D steady incompressible flow inside a triangu- 
lar cavity is governed by the Navier-Stokes equations. We 
consider the N-S equations in streamfunction and vorticity 
formulation, such that 
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Substituting for the transformation metrics defined in 
Equations (3 and 4), we obtain the equations that govern 
the flow in a triangular cavity shown in Figure 1 as the 
following 
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where 
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The boundary conditions in x- and j/-plane arc the 
following 

ip = on three sides 
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on the top side 
on the other two sides 

on three sides 
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where n is the unit normal vector and t is the tangent 
vector in clockwise direction and (u, v) are the velocity 
components in x- and y-direction where 
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Therefore, in the computational domain, the boundary 
conditions on the side A-B are 



dip = dip = dy 

d£ drj drj 

On sides A-C and B-C we have 



dip 
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Also note that, on side B-C one can also show that the 
derivative of the streamfunction in the wall normal direc- 
tion would be homogenous. 

Erturk et al. [6] stated that for square driven cavity 
when fine grids are used, it is possible to obtain numerical 
solutions at high Reynolds numbers. At high Reynolds 
numbers thin boundary layers are developed along the 
solid walls and it becomes essential to use fine grid meshes. 
Also when fine grids are used, the cell Reynolds number 
or so called the Peclet number, defined as Re c = uAh/v, 
decreases and this improves the numerical stability (see 
Weinan and Jian-Guo [22] and Tennehill et al. [20]). In 
this study, we used a fine grid mesh with (512x512)/2 
grid points. On this mesh, we solved the governing equa- 
tions (10 and 11) using SOR method. For the vorticity 
values at the boundary points we used Thorn's formula 
[21]. We note an important fact that, it is well under- 
stood (see Weinan and Jian-Guo [22], Spotz [19], Huang 
and Wetton [10], Napolitano et al. [16]) that, even though 
Thorn's method is locally first order accurate, the global 
solution obtained using Thorn's method preserve second 
order accuracy 0{Ah 2 ). On side A-B the vorticity is equal 
to 
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where the subscripts i and j are the grid indices in £- 
and 77-directions respectively and the grid index refers 
to points on the wall and 1 refers to points adjacent to the 
wall and also Ah refers to grid spacing. 

Similarly on side A-C the vorticity is equal to 
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and on side B-C, it is equal to 
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3 RESULTS AND DISCUSSIONS 



During the iterations, as a measure of convergence, we 
monitored several residuals. We define the residual of the 
steady state streamfunction and vorticity Equations (10 
and 11), as the following 
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These residuals indicate the degree to which the nu- 
merical solution has converged to steady state. In our 
computations, in all of the triangle geometries considered, 
at any Reynolds numbers, we considered that the conver- 
gence was achieved when the absolute values of the maxi- 
mum of the residuals defined above in the computational 
domain (max(|i?^|) and max(|i? w |)) were both below 
10 -10 . Such a low convergence level is chosen to ensure 
the accuracy of the solution. At these convergence levels, 
the maximum absolute change in the streamfunction vari- 
able in the computational domain (max( lip™^ 1 — ip™j \ )) 
was in the order of 10 -16 and for the vorticity variable 
(maxda^t 1 — it was in the order of 10~ 14 . In 

other words, our numerical solutions of the streamfunc- 
tion and the vorticity variables are accurate up to 15 and 
13 digits respectively. Also at these convergence levels, the 
maximum absolute normalized change in the streamfunc- 
tion and vorticity variable in the computational domain 
(max( m^-^W/^U ) andmax( K^t 1 -^)/^! )) 
were in the order of 10~ 13 and 10~ 12 respectively. In other 
words, at convergence the streamfunction variable was 
changing with 10 -11 percent of its value and the vorticity 
variable was changing 10~ 10 percent of its value, at each 
iteration step. These numbers ensure that our numerical 
solutions are indeed very accurate. 
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In this study we considered different triangle geome- 
tries. In order to demonstrate how we define the Reynolds 
number, let us consider a dimensional equilateral triangle 
with coordinates of corner points 

a x = —VSa , b x = V3a , h — a , c x = , c y = —2a (21) 

If we use the top wall velocity U and the length scale 
a in nondimcnsionalization, then the Reynolds number is 
defined as Ua/v and the coordinates of the corner points 
of the nondimensional equilateral triangle become 

a x = -V3 , b x = V3 , h = 1 , c x = , c y = -2 (22) 

We, first considered this equilateral triangle geometry 
which was also considered by McQuain et al. [14], Ribbens 
et al. [17], Jyotsna and Vanka [11], Li and Tang [13] and 
Gaskell et al. [8]. We solved the flow in this triangle cav- 
ity at various Reynolds numbers ranging between 1 and 
1750. We note that, if the length of one side of the trian- 
gle (2>/3fl) was used in nondimcnsionalization, as it was 
used by Gaskell et al. [8], then our Reynolds number of 
1750 would be 2\/3 fold such that it would correspond to 
a Reynolds number of 6062. Qualitatively, Figure 2 shows 
the streamline contours of the flow at various Reynolds 
numbers and also Figure 3 shows the vorticity contours 
at the same Reynolds numbers. In terms of quantitative 
analysis, Table 1 tabulates the center of the primary eddy 
and the streamfunction and vorticity values at the core, 
together with results found in the literature. Our results 
agree well with McQuain et al. [14], Ribbens et al. [17] 
and agree excellent with Gaskell et al. [8] up to the max- 
imum Reynolds number (i?e=500) they have considered. 
However, the results of Li and Tang [13] differ from our 
results. The difference can be seen more obvious in Fig- 
ure 4, where we plot the vorticity values at the center 
of the primary vortex tabulated in Table 1, with respect 
to the Reynolds number. The results of Li and Tang [13] 
start to behave differently starting from i?e=100 from the 
rest of the results. Also in Figure 4, the results of Mc- 
Quain et al. [14] and Ribbens et al. [17] start to devi- 
ate from present results and Gaskell et o/.'s [8] results 
after i?e=200. Having a different behavior, the vorticity 
value of McQuain et al. [14] and Ribbens et al. [17] shows 
an increase such that their vorticity value at i?e=500 is 
greater than the vorticity value at i?e=350, whereas the 
present results and Gaskell et al.'s [8] results show a con- 
tinuous decrease. We believe that these behaviors are due 
to the coarse grids used in those studies and in order to 
resolve these behaviors we decided to solve the same flow 
problem using several coarse grid meshes. We computed 
the flow using (64x64)/2, (128xl28)/2 and (256x256)/2 
grid points also. For a given grid mesh we have solved the 
flow for increasing Reynolds number until a particular Re 
where the solution was not convergent but oscillating. For 
this particular Re when the number of grid points was in- 
creased, the convergence was recovered and we were able 
to obtain a solution. 



In order to demonstrate the effect of the grid mesh 
on the numerical solution, in Figure 5 we plotted the 
vorticity value at the center of the primary vortex ob- 
tained using different grid mesh, with respect to Reynolds 
number. From the results in Figure 5, we conclude that 
(64x64)/2 grid mesh is too coarse for this flow problem. 
In this figure, looking at the solution of (128xl28)/2 and 
(256x256)/2, the vorticity value decreases with increas- 
ing Reynolds number however shows an increase at the 
largest Reynolds number where beyond that Re the so- 
lution is not converging. We note that this behavior is 
similar to that of McQuain et al. [14] and Ribbens et al. 
[17] in Figure 4. 

In Figure 5 looking at the solutions of different grid 
mesh, for example at i?e=350, we see that as the number 
of grid points in the mesh is increased, the vorticity value 
at the center of the primary eddy is decreased however this 
decrease slows down as Ah gets smaller. This is not sur- 
prising since the solution gets more accurate as the num- 
ber of grid points is increased, and also since our solution 
is spatially second order accurate, 0(Ah 2 ), the computed 
vorticity value should be expected to converge to the real 
value with Ah 2 . At the finest grid mesh, (512x512)/2, at 
high Reynolds numbers our the vorticity value seems to 
stagnate around 1.16. We note that, we expect the vortic- 
ity value to be a little lower than this value, if the number 
of grid points is increased furthermore. However increas- 
ing the number of grid number furthermore increases the 
computing time, therefore we will not consider a finer grid 
mesh than (512x512)/2. This grid study shows that fine 
grids are necessary for accurate solutions for triangular 
cavity flow especially at high Reynolds numbers. 

McQuain et al. [14] analytically obtained the constant 
vorticity value of the inviscid rotational flow at infinite 
Reynolds number as vT0/3=1.054 using Batchelor's mean- 
square law. This infinite Reynolds number value of the 
vorticity of 1.054 is shown with the dotted line in Fig- 
ure 4. Our computed primary eddy vorticity asymptotes 
to a value around 1.16 as the Reynolds number increases, 
which is a little higher than the theoretical value. We note 
that the theory, assumes that at infinite Reynolds num- 
ber the whole inviscid fluid rotates as a solid body with 
a constant vorticity which is coupled to boundary layer 
flows at the solid surface. However in reality, from Fig- 
ure 2 it is evident that the whole fluid does not rotate as 
a solid body, instead there appears progressively smaller 
counter-rotating recirculating regions at the bottom cor- 
ner and also one towards the upper left corner. We note 
that, especially the eddies at the bottom corner occupy 
some portion of the corner as the Reynolds number in- 
creases. However, the increase in the size of the portion of 
the bottom corner eddies almost stops after i?e=500 and 
the size of the primary eddy remains almost constant be- 
yond and so is the value of the vorticity at the core of the 
primary eddy. It looks like, for an equilateral triangular 
cavity flow at high Reynolds numbers, the mean square 
law predicts the strength of the primary eddy within an 
error due to the effect of the secondary eddies. For circular 
or elliptic boundaries Ribbens et al. [17], for square cavity 
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flow Erturk et al. [6], and for rectangle cavities McQuain 
et al. [14] have shown that the mean square law is approx- 
imately valid and successful in predicting the strength of 
the primary eddy at high Reynolds numbers. Due to small 
corner angles in a triangle geometry, for triangular cavi- 
ties the mean square law is not as successful as it was in 
other geometries. 

We then considered an isosceles triangle which was also 
considered by Jyotsna and Vanka [11], where 

a x = -1 , b x = 1 , h = , c x = , c v = -4 (23) 

We note that, our definition of Reynolds number is 
equivalent to one fourth of the Reynolds number definition 
used by Jyotsna and Vanka [11]. 

Figure 6 shows the streamline contours of the flow in 
this triangle at various Reynolds numbers. Table 2 com- 
pares the location of the primary eddy center with that of 
Jyotsna and Vanka [11] and Gaskcll et al. [8]. The results 
agree with each other, although we believe that our results 
are more accurate. 

Moffatt [15] analytically studied the recirculating ed- 
dies near a sharp corner in the Stokes regime and predicted 
that the distance from the corner to the centers of eddies 
and also the velocity at the dividing streamline between 
the eddies, follow a geometric sequence. 

Using the considered triangle geometry, Jyotsna and 
Vanka [11] have used the small Reynolds number solutions 
to compare with Moffatt's [15] Stokes regime predictions. 
However, when Reynolds number increases the asymme- 
try in the flow field increases and the comparison with 
the Stokes regime could not be possible. In order to be 
able to compare our results with Moffatt [15], we solve 
the considered triangle geometry for the Stokes regime, 
i.e. Re=0. The solution at Re=0, as seen in Figure 6a, is 
symmetrical. For this Stokes regime solution, we first cal- 
culate the w-velocity along a perpendicular line from the 
bottom corner to the top wall. In this velocity profile, the 
points where the u-velocity is zero correspond to the eddy 
centers. The distance of these points to the bottom corner 
would give us r n , where we use in calculating the ratio 
fn/fn+i- I n the velocity profile, the points where the it- 
velocity has local maximum correspond to the ^-velocity 
at the dividing streamline between the eddies, where Mof- 
fatt [15] have used these velocities as a measure of the 
intensity of consecutive eddies. These local maximum u- 
velocities would give us /„, where we use in calculating 
the ratio I n /I n+ \. For the considered triangle geometry 
where 0=28.072°, Table 3 tabulates the calculated ratios 
of r n /r n+ i and I n /I n +i along with analytical predictions 
of Moffatt [15] and the agreement is good. 

We then considered an isosceles triangle which was also 
considered by Gaskell et al. [8], with 

a x = -1 , b x = 1 , h = , c x = , c y = (24) 

tan(|) 

where (9=20°, 30°, 40°, 50°, 60°, 70°, 80° and 90°. We will 
consider Stokes flow in these triangles, such that we solve 



for Re=0. The analytical predictions of Moffatt [15] for 
these 9 values, arc tabulated in Table 4. Figure 7 qualita- 
tively shows the change in the Stokes flow in an isosceles 
triangle as the corner angle, 9, changes. For this case, our 
computed relative positions and intensities of eddies are 
tabulated in Table 5 along with results of Gaskell et al. 
[8] . The agreement between our computed results and that 
of Gaskell et al. [8] and also predictions of Moffatt [15] is 
good. We note that, Moffatt's [15] analytical analysis pre- 
dicts that there will be infinite sequence of eddies at the 
bottom corner. This result is, of course, impossible to cap- 
ture for a numerical analysis due to the resolution of the 
grid mesh. However since we have used a finer grid mesh, 
with (512x512)/2 grid points, than the grid mesh used by 
Gaskell et al. [8], with 6605 nodes, we were able to cap- 
ture more sequence of eddies than Gaskell et al. [8] , at the 
sharp bottom stagnant corner. 

We then considered an isosceles right triangle with the 
90° corner being at the top right corner, such that with 
corner points 

a x =0,b x = l,h = l, c x = l , c y = (25) 

The flow in this considered triangle is solved up to 
i?e=7500. Figure 8 shows the flow topology as the Reynolds 
number increases and Table 6 tabulates the properties of 
the primary eddy, the streamfunction and the vorticity 
value at the center of the eddy and also the location of 
the center, for future references. 

We then also considered an isosceles right triangle with 
the 90° corner being at the top left corner, such that with 
corner points 

a x = , b x = 1 , h = 1 , c x = , c y = (26) 

Using this triangle geometry, we were able to obtain 
solutions up to i?e=2500. Figure 9 shows the flow topol- 
ogy as a function of Reynolds numbers. In Table 7, for 
the considered triangle geometry, we tabulated the prop- 
erties of the primary eddy, i.e. the eddy that is closest to 
the moving lid, for future references. As it is obvious from 
Figures 8 and 9, in both cases the flow behaves very differ- 
ently as the Reynolds number increases, which shows that 
the flow structures in a triangle cavity are greatly affected 
by the triangle geometry. 

4 CONCLUSIONS 

In this study we have presented high accurate, fine grid 
solutions of 2-D steady incompressible flow in triangle cav- 
ities. The governing equations are solved up to very low 
residuals at various Reynolds numbers. The computed so- 
lutions are compared with numerical solutions of McQuain 
et al. [14], Ribbcns et al. [17], Jyotsna and Vanka [11], Li 
and Tang [13] and Gaskcll et al. [8] and also with ana- 
lytical predictions of Moffatt [15], and good agreement is 
found. Our results showed that for an equilateral trian- 
gular cavity flow, Batchelor's mean-square law is not as 
successful as it was in square or rectangle cavity flows, 
due to small stagnant corner angle. 
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Figure 1) Geometric transformation of the triangular cavity 
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Figure 2) Streamline contours at various Reynolds numbers 
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Figure 3) Vorticity contours at various Reynolds numbers 




Figure 4) Comparison of the vorticity at the center of the primary eddy 




Figure 5) Vorticity at the center of the primary eddy obtained with different grid mesh 




d) Re=100 e) Re=200 



Figure 6) Contour figures of isosceles triangle 




Figure 8) Contour figures of right 



hand side aligned right triangle at various Re. 




Figure 9) Contour figures of left hand side aligned right triangle at various Re. 
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1 
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Table 1) Comparison of primary eddy properties, for equilateral triangle 
with corner points a x = -y[3 , b x = -J3 , h = 1 , c x =0, c = -2 
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(0.225,-0.450) 
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Table 2) Comparison of the location of the primary eddy center for isosceles 
triangle with corner points a x = -1 , b x =l, h = 0, c y = , c y = -4 
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Table 3) Relative eddy center locations r n I r n+1 and intensities I n II. 
for isosceles triangle with 9 = 28.072° 



e 




I J In* 


20° 


1.64 


381.5 


30° 


2.11 


414.5 


40° 


2.76 


467.9 


50° 


3.67 


552.2 


60° 


5.00 


687.4 


70° 


7.05 


915.2 


80° 


10.44 


1331.1 


90° 


16.57 


2189.1 



Table 4) Moffatt's [15] predictions of relative eddy center locations r n I r t 
and intensities I n I I n+1 for Stokes regime for various 6 
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Re 
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(x,y) 


100 


-0.07162 


-4.83000 
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Table 6) The properties of the primary eddy, for right hand side aligned right triangle 
as a function of Reynolds number 



Re 
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CO 


(x,y) 
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-0.06451 
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-0.05306 


-7.02235 
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Table 7) The properties of the primary eddy, for left hand side aligned right triangle 
as a function of Reynolds number 



